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Reliable calculations of nuclear binding energies are crucial for advancing the research of nuclear physics. 
Machine learning provides an innovative approach to exploring complex physical problems. In this study, the 
nuclear binding energies are modeled directly using a machine-learning method called the Gaussian process. 
First, the binding energies for 2238 nuclei with Z > 20 and N > 20 are calculated using the Gaussian process 
in a physically motivated feature space, yielding an average deviation of 0.046 MeV and a standard deviation of 
0.066 MeV. The results show the good learning ability of the Gaussian process in the studies of binding energies. 
Then, the predictive power of the Gaussian process is studied by calculating the binding energies for 108 nuclei 
newly included in AME2020. The theoretical results are in good agreement with the experimental data, reflect- 
ing the good predictive power of the Gaussian process. Moreover, the a-decay energies for 1169 nuclei with 
50 < Z < 110 are derived from the theoretical binding energies calculated using the Gaussian process. The 
average deviation and the standard deviation are, respectively, 0.047 MeV and 0.070 MeV. Noticeably, the cal- 
culated a-decay energies for the two new isotopes 7°* Ac [M. H. Huang et al., Phys. Lett. B 834, 137484 (2022)] 
and ??"Th ІН. B. Yang et al., Phys. Rev. C 105, L051302 (2022)] agree well with the latest experimental data. 
These results demonstrate that the Gaussian process is reliable for the calculations of nuclear binding energies. 
Finally, the a-decay properties of some unknown actinide nuclei are predicted using the Gaussian process. The 
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predicted results can be useful guides for future research on binding energies and a-decay properties. 
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I. INTRODUCTION 


Nuclear binding energies are important ground state prop- 
erties that provide valuable information for probing nuclear 
structures [1—4] and serve as crucial inputs for some nuclear 
physics problems [5, 6]. For instance, binding energies play 
a key role in calculating the product cross sections for un- 
known nuclei using nuclear reaction models before synthe- 
sizing superheavy nuclei [7, 8]. They are also instrumen- 
tal in identifying new nuclides in synthesis experiments of 
heavy and superheavy nuclei [9, 10] because o decay is one 
of the fundamental decay modes for most heavy and super- 
heavy nuclei [11—13]. For o-emitters, there are two main a- 
decay observable properties, which are respectively a-decay 
energies and half-lives [14—18]. Thereinto, a-decay half-lives 
are strongly influenced by the a-decay energies, which can 
be calculated using the binding energies. Meanwhile, bind- 
ing energies are also vital for calculating the properties of 
other radioactive decay modes, such as two-proton radioac- 
tivity [19] and heavy-cluster radioactivity [20]. Furthermore, 
the accuracy of binding energies has a significant impact on 
nuclear astrophysics studies, including r-process [21, 22], rp- 
process [23, 24], and the properties of neutron stars [25, 26]. 
Therefore, it is necessary to explore reliable theoretical mod- 
els to calculate and predict the binding energies more accu- 
25 rately. 
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With the advancements in experimental nuclear physics fa- 
cilities, binding energies of more than two thousand nuclei 
have been measured to date [27]. The accumulated experi- 
mental data provide a foundation for the development of the- 
oretical models. In the past few years, numerous theoret- 
ical models and formulas have been proposed to calculate 
binding energies, including the Bethe-Weizsacker formula 
[28, 29], the Thomas-Fermi (TF) model [30], the Hartree- 
Fock-Bogoliubov mean field model [31], and the finite-range 
drop model (FRDM) [32]. The theoretical binding energies 
calculated using these models and formulas are in good agree- 
ment with the experimental data. In Ref. [8], an improved 
binding-energy formula was proposed by incorporating ad- 
ditional physical terms into the standard Bethe-Weizsacker 
formula, which consists of the shell effect and the neutron- 
proton correlations. The binding energies and a-decay ener- 
gies can be well reproduced using this improved formula for 
heavy and superheavy nuclei with Z > 90 and N > 140. 
Although these current traditional models can provide theo- 
retical guidance for studying binding energies, it is still worth 
exploring other models to provide more accurate calculations 
and predictions for future investigations of binding energies. 


Machine learning has been widely used across many fields 
[33—38], as it can learn useful information from known sys- 
tems and predict unknown properties within the same system 
using the obtained information. In the last decade, nuclear 
properties have been studied using various machine-learning 
methods based on available physical knowledge, including 
nuclear masses [39—41], nuclear charge radii [42], a-decay 
properties [43], and 8-decay properties [44]. These nuclear 
properties can be well reproduced using machine learning. 
Recently, a new Bayesian machine learning mass model has 
been proposed [45], which can reproduce nuclear masses with 
the high accuracy required for the studies of r-process. As 
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one of the popular machine-learning methods, the Gaussian 
process is a powerful nonparametric model, which is expected 
to model any distribution of the objectives [46]. Owing to 
its excellent flexibility in data modeling, the Gaussian pro- 
cess has been frequently applied in various studies [47, 48]. 
Notably, the Gaussian process can provide not only the the- 
oretical values of the objectives but also the distribution of 
the calculated results, contributing to the visualization of the 
theoretical uncertainties [49]. Recently, the Gaussian process 
has been successfully exploited to predict the a-decay ener- 
gies and half-lives of actinide nuclei [50]. Inspired by these 
previous works, it is of great interest to explore the reliability 
of the Gaussian process in the calculations of binding ener- 
gies. 

In this work, the Gaussian process has been extended to 
study the binding energies by directly modeling the experi- 
mental binding energies. The remainder of this paper is given 
as follows. In Sec. II, the theoretical framework, consisting 
of the Gaussian process with the modified kernel function 
and the physically motivated feature space, is provided. In 
Sec. III, the theoretical binding energies calculated using the 
Gaussian process are shown and discussed. Furthermore, the 
a-decay properties are reproduced and predicted based on the 
calculated binding energies. Finally, a comprehensive sum- 
mary is presented in Sec. IV. 


П. THEORETICAL FRAMEWORK 


In the present work, the binding energy for a nucleus is 
considered as a realistic observation By = bp + 6 with 
noise 0 ~ (0,02). Here, bj = b(x,) is a latent func- 
tion that denotes the noise-free binding energy for the pth 
nucleus x, [51]. В, denotes the realistic binding energy, 
and ð is an independently identically distributed Gaussian 
noise. Given a set of n nuclei with known binding ener- 
gies into a training set (xp, В)! 1, we aim to model the 
underlying physical relationship between each nucleus and 
its binding energy using the Gaussian process. Within the 
framework of the Gaussian process, the values of latent func- 
tion b = (b1, ba, +++ , bn)? = (b(x1), b(x2), +++ bis, IT are 
modeled by a joint Gaussian distribution, characterized by 
the values of a mean function (m(x1), m(x2),--- , M(Xn))T 
and the matrix of a covariance function [k(x,, Xq)]nxn [46]. 
Therefore, the Gaussian process can be generally denoted as 
(хы) ~ GP(m(x,), k(x,, x,)). The mean function m(x,) 
is often set as zero because of the lack of prior knowledge. 
The so-called kernel function k(xp, x4) can be written as а 
function of |x, — Sol, which is crucial for describing the sim- 
ilarities between pairs of nuclei. For the studies of binding 
energies, we choose a composite kernel function written as 


k(xy, X4) = т 


with ть = |х, — x4|. The modified kernel function is a linear 
combination of two widely used kernel functions, which are 
the Matérn kernel function and the Rational Quadratic ker- 
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nel function, respectively. Here, ть, ly, оь, and dy are four 
hyperparameters of the Gaussian process. lp, оь, and dr can 
capture the relevant range of the binding energies for pairs 
of nuclei, and m is able to describe the correlation intensity 
between them. For the realistic binding energy Bp, the co- 
variance function becomes k(xy, хи) > k(Xp,Xq) + о?б,а. 
бра is a Kronecker delta where Au = 1 for p q and 
Au = 0 for p # а. When describing a number of nu- 
clei X = (xi,Xo,--,x,)7, the binding energies B 
(B4, Bo,+-+ , Bn)" are expressed as B ~ GP(0, K(X, X)+ 
cL), where I is a diagonal matrix. 

The central interest of this work is to predict unknown 
binding energies based on the knowledge learned from the 
training set using the Gaussian process. When predicting un- 
known binding energies for nuclei X. with the training set 
D = (Xp, By) 1, the joint Gaussian distribution of the train- 
ing outputs B and the predicted outputs b, can be written as 


[46] 
B 

al ) е 
For n, predicted nuclei, K(X, X), K(X, X.), K(X., X), 
and K(X., X..), respectively, denote n x n, n x Nx, n. x n, 
and n, x n, matrix evaluated at all pairs of training and pre- 
dicted points. By conditioning the joint Gaussian distribution, 
the crucial predicted expressions for the Gaussian process are 


b.[D, X. ~N (bx, cov(b.)), where 


K(X,X)--o?I K(X, X.) 
K(X,X) қ(х,,х.) 


b, - K(X,, X)[K(X, X) + o?T|-! B, 


cov(b.) =K(X,, X,) - K(X,, X)[K(X, X)-o?I] ! K(X, X.). 


(3) 
Here, the values of b, give the predicted binding energies 
for unknown nuclei. The variances of the predicted binding 
energies can be calculated by adding the noise variance о? to 
the predictive variance given by соу(Ь,). 

As mentioned above, each nucleus is described by x,, 
which is a vector of physical features determining the de- 
scription of the corresponding binding energy. In the present 
work, our goal is to obtain good descriptions of the bind- 
ing energies using the Gaussian process with as simple 
physical information as possible. Hence, we construct a 
physical feature space with nine features, where x, = 
(Ap, AJ", gae (Ap/2 ЕЕ 2,)2/4;, Ag Ops ДА SÉ 
Zy|/ Ap, Tp, Up). Неге, A, Z, and N denote the mass, proton, 
and neutron numbers, respectively. The first six features are 
based on the Bethe-Weizsacker formula [7, 28, 29, 52, 53]. A 
is introduced to model the proportional relationship between 
the binding energies and the nuclear volume, reflecting the 
saturation of nuclear force. A?/* is provided since the bind- 
ing energies are expected to decrease on the nuclear surface. 
Z2 A-1/? is used to describe the influence of the Coulomb 
interaction between protons. (4/2 — Z)?/A is the symmetry 
term that approximately estimates the balance between N and 
Z. A-1? and 6 = [(—1)N + (—1)7]/2 are used to describe 
the pairing energies with ó — 1,0, —1 for the even-even, odd- 
A, and odd-odd nuclei. |N — Z|/A is from the Wigner term, 
which originates from the neutron-proton correlations [1, 8]. 
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Additionally, т and v include the shell information, where 
7 (v) is calculated using the numbers of protons (neutrons) 
away from the nearest proton (neutron) magic numbers [54]. 

The aforementioned theoretical framework implies that 
five hyperparameters need to be determined, which are 7p, lp, 
оь, dp, and a», respectively. These сап be determined by op- 
timizing the marginal likelihood using the training data [46]. 


Ш. NUMERICAL RESULTS AND DISCUSSIONS 


In this section, we present and discuss the theoretical re- 
sults of the nuclear binding energies calculated using the 
Gaussian process. First, we calculate the binding ener- 
gies for nuclei with Z > 20 and N > 20 to evaluate 
the learning ability of the Gaussian process. The training 
set chosen in this work contains 2238 nuclei with known 
binding energies taken from AME2020 [55]. Each nucleus 
in the training set is presented as (xp, Bp), where x, = 
(Ap, АШ АСТА — 2 DA, Ap’, õp, |Np — 
Zp|/Aps Tp, Vp) and B, = BE*P*. After the training pro- 
cess, the hyperparameters are determined as 7), = 1.814 x 
104 MeV!/2, lp = 1.821 x 104, ap = 1937.218, dy = 
414.771, and оу = 0.093 MeV 2, The larger value of rj, in- 
dicates a stronger dependence between pairs of nuclei. Mean- 
while, the larger values of l, оь, and d; result іп a relatively 
larger correlation range, which means that the change of bind- 
ing energies is comparatively smoother. Moreover, they also 
assist in avoiding the rapid growth of the error bars of the 
binding energies for nuclei away from the training data [46]. 

After the hyperparameters have been determined, the bind- 
ing energies can be calculated using the Gaussian process. To 
test the accuracy of the calculated results, we calculate the ab- 
solute value of the deviation between the experimental result 
and the theoretical one for each nucleus, defined by 

ІА в| — Ipse _ Т . (4) 
Here, Bre and ВІ һео- denote the experimental binding еп- 
ergy and theoretical result calculated using the Gaussian pro- 
cess for the pth nucleus, respectively. The numerical results 
show that all absolute values of the deviations are smaller than 
0.423 MeV, indicating a small global deviation. We show the 
corresponding results in Fig. 1, in which the x- and y-axis in- 
dicate the neutron and proton numbers, respectively. The red 
squares depict the absolute values of the deviations, where 
darker colors are associated with larger deviations. The trans- 
verse and vertical dotted lines present № = 28, 50, 82, 126 
and Z = 28, 50,82, respectively. It can be seen clearly from 
Fig. 1 that the colors of most squares are lighter, reflecting 
that the deviations for most nuclei are below 0.1 MeV. Addi- 
tionally, the binding energies for nuclei near the shell closure 
are also well reproduced using the Gaussian process. Next, 
we calculate the average deviation 
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Fig. 1. Тһе absolute values of deviations between experimen- 
tal binding energies and the theoretical results calculated using the 
Gaussian process across the nuclear chart. The darker colors indi- 
cate larger deviations of binding energies. Numerically, the largest 
absolute value of the deviations is |А в| = 0.423 MeV. 


and the standard deviation 
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of the theoretical binding energies calculated using the Gaus- 
sian process for nuclei with Z > 20 and N > 20. Here, 
Ав denotes the number of nuclei included in the calcula- 
tions. The numerical values аге (св) 0.046 MeV and 
oy 0.066 MeV, respectively. The small deviations 
show that the theoretical binding energies calculated using 
the Gaussian process with the modified kernel function in 
the physically motivated feature space are in good agreement 
with the experimental data. These results demonstrate the 
good learning ability of the Gaussian process in the studies 
of binding energies. 

To further evaluate the learning ability and predictive 
power of the Gaussian process in the studies of binding en- 
ergies, we perform cross validation for the Gaussian process. 
In this work, we introduce the isotone-fold cross-validation 
that nuclei in each isotonic chain will be predicted using the 
Gaussian process based on the information provided by the 
remaining isotonic chains in the training set. The average de- 
viations and the standard deviations of the theoretical binding 
energies for nuclei in each isotonic chain are calculated, with 
results depicted in Fig. 2. For comparison, the average devia- 
tions and the standard deviations of the binding energies cal- 
culated using the Bethe-Weizsacker formula for each isotonic 
chain are also provided in Fig. 2. In Fig. 2(a) and Fig. 2(b), 
the red squares denote the average deviations and the standard 
deviations calculated using the Gaussian process for each iso- 
tonic chain, respectively. The blue circles present the average 
deviations and the standard deviations calculated using the 
Bethe-Weizsacker formula for each isotonic chain separately. 
It is straightforward to see that the deviations given by the 
Gaussian process are quite small, which means that the cross- 
validation result is pretty good. In addition, we can find that 
the deviations are significantly reduced compared with those 
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The cross-validation results for nuclei in each isotopic chain calculated using the Gaussian process. In Fig. 2(a), the red squares 


and the blue circles depict the average deviations calculated using the Gaussian process and the Bethe-Weizsacker formula, respectively. In 
Fig. 2(b), the red squares and the blue circles show the standard deviations calculated using the Gaussian process and the Bethe-Weizsacker 


formula separately. 


given by the Bethe-Weizsacker formula. These results reflect 
the good learning ability and predictive power of the Gaussian 
process. Numerically, the total average deviation and stan- 
dard deviation of the cross-validation for nuclei in the train- 
ing set are (op) 0.100 MeV and 4/02 = 0.144 MeV, 
respectively. The small deviations show that the predicted 
binding energies agree well with the experimental data, in- 
dicating that the binding energies can be well learned using 
the Gaussian process. Thus, we can conclude that the learn- 
ing ability and predictive power of the Gaussian process are 
reliable for studying the binding energies. 

Then, we further test the predictive power of the Gaussian 
process by calculating the binding energies for nuclei that are 
present in AME2020 but not in AME2012 using the Gaussian 
process. To perform this calculation, the training set is cho- 
sen to include nuclei that are provided in both AME2012 and 
АМЕ2020. Based on the training set, we predict the bind- 
ing energies for 108 nuclei that are provided in AME2020 but 
not in AME2012 using the Gaussian process. The theoretical 
average deviation and standard deviation for these nuclei are 
(cg) = 0.216 MeV and Vo 0.304 MeV, respec- 
tively. These deviations are acceptable results in the calcula- 
tions of binding energies, verifying that the predicted power 
of the Gaussian process is commendable. Therefore, based on 
these theoretical results, it can be concluded that the Gaussian 
process is a reliable model for the studies of nuclear binding 
energies. 

Next, we would like to calculate and discuss the theoret- 
ical results calculated using the Gaussian process with dif- 
ferent kernel functions and physical feature spaces. First, 
we calculate the binding energies using the Gaussian pro- 
ғо cess with the Matérn kernel function and the Rational 
231 Quadratic kernel function, respectively. Тһе correspond- 
zs ing deviations of the binding energies for 2238 nuclei are 
әз (oni, / 02) (0.059, 0.076) MeV for the Matérn ker- 
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әм nel function and ((ов), \/o%) (0.121, 0.166) MeV 
285 for the Rational Quadratic kernel function, respectively. 
ae The deviations for 108 new nuclei аге ((op) , у/о2) 
27 (0.278, 0.415) MeV for ће Matérn kernel function and 
oe (op). / 025) (0.193,0.249) MeV for the Rational 
239 Quadratic kernel function separately. Comparing with the de- 
о Viations ((5) , (/02,) = (0.046,0.066) MeV for 2238 nu- 


ә clei and ((св), vei) = (0.216, 0.304) MeV for 108 new 
2»» nuclei calculated using the composite kernel function, it can 
23 be found that the deviations calculated with the composite 
294 kernel function are as small as those calculated using the 
295 Matérn kernel function for 2238 nuclei and show better re- 
296 Sults than those calculated using the Matérn kernel function 
297 for 108 new nuclei. The deviations calculated using the com- 
ep posite kernel function show results as good as those calcu- 
lated using the Rational Quadratic kernel function for 108 
new nuclei and are smaller than those calculated using the Ra- 
tional Quadratic kernel function for 2238 nuclei. Therefore, 
the good interpolation power of the Gaussian process with 
the Matérn kernel function and extrapolation ability of the 
Gaussian process with the Rational Quadratic kernel function 
are inherited by the composite kernel function in the calcula- 
tions of binding energies, which demonstrates that the mod- 
ified kernel function is a good choice for the present work. 
зов Furthermore, we hope that the choice of the composite kernel 
зө function can provide a new idea for modeling other physical 
зо problems using the Gaussian process. 
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We continue to compare the average deviations and stan- 
dard deviations for 108 new nuclei using the Gaussian pro- 
cess in different physical feature spaces. We first cal- 
culate the deviations for nuclei using the Gaussian pro- 
cess in the feature space consisting of six features taken 
from the Bethe-Weizsacker formula, where the pth nu- 


cleus is described by x, = (Ay, A2”, ZA”, (4,/2- 


318 


3 


9 


32 


o 


32 


322 


323 


324 


32! 


a 


32 


o 


32 


N 


32i 


© 


32! 


© 


33l 


б 


33 


332 
333 
334 


33! 


a 


336 


33 


e 


33i 


Ф 


33! 


Ф 


34 


о 


34 


342 
343 
344 


ЕЛ 


a 


346 


34 


E 


348 
34! 


© 


35 


о 


35 


352 
35 


[5] 


354 


Zo) As Ap e 5). The theoretical deviations аге 
(oni, vol = (0.437,0.775) MeV. Then, we add the 
neutron-proton correlation and the shell information in the 
above feature space and compare the corresponding devia- 
tions. When the neutron-proton correlation is added in the 
feature space where x, — (LE gea ue (A,/2 — 
Zp)? /Ap, Ap”, õp, |Np — Zp|/Ap), the deviations become 
(ов), \/o%) = (0.398, 0.712) MeV. The reduction in the 
deviations shows that the neutron-proton correlation is nec- 
essary for calculating the binding energies. When the shell 
information is included in the feature space, where x, = 
2/3 2 4-1/3 2 -1/2 

(Ар, Ар , Zp Ар , (Ap/2 — Zp) /Ар, Ар , Ôp, Tp, Up), 
the deviations аге (св) , y/o} ) = (0.236, 0.365) MeV. The 
results reflect that the introduced features 7 and v provide 
useful shell information for nuclei in the calculations of bind- 
ing energies. Furthermore, it can be observed that the above 
deviations are larger than those calculated in the feature space 
with nine features established in the present work, indicating 
that our choice of feature space is reasonable. Notably, the 
importance of the physically motivated feature space has also 
been studied in the Bayesian neural network and the proba- 
bilistic Mixture Density Network [39, 41]. The physical fea- 
ture space established in the present work is first studied in 
the Gaussian process on the research of binding energies. 

It has been mentioned that the distribution of theoretical 
results can be provided by the Gaussian process. Here, we 
present the intervals of error bars for the theoretical results 
calculated in this work. The lengths of error bars at 95% 
confidence interval range from 0.213 MeV to 0.258 MeV in 
the studies of 2238 nuclei, while they range from 0.234 MeV 
to 4.022 MeV in the calculations of 108 new nuclei. These 
results show that the hyperparameters determined by the 
marginal likelihood are reasonable and that the theoretical 
binding energies calculated using the Gaussian process are 
reliable. Thus, we conclude that the Gaussian process with a 
modified kernel function and the physically motivated feature 
space is a reliable model for calculating binding energies. 


Table 1: The theoretical a-decay energies calculated using the Gaus- 
sian process for some actinide nuclei. The first column denotes the 
actinide nuclei. The second and third columns list the experimen- 
tal o-decay energies and the theoretical values calculated using the 
Gaussian process separately. The last column presents the deviations 
AQ, = ОЕ" — QTh°°-. The experimental data for the new nuclides 
204 Ac and 7°" Th are taken from Ref. [57] and Ref. [10], respectively. 


Nucl. QEPE (MeV) Theo. (MeV) AQa (MeV) 
204 Ac [57] 8.107 8.107 0.000 

205 Ac 8.093 8.083 0.010 

206 Ac 7.958 7.943 0.015 

207 Ac 7.845 7.863 -0.018 

208 Ac 7.729 7.136 -0.007 

209 Ac 7.730 7.103 0.027 

210 A c 7.586 7.608 -0.022 

211 Ac 7.568 7.569 -0.001 


Table 1: (continued) 


Nucl. QEPE: (MeV) Theo. (MeV) AQ. (MeV) 
ZIZ Ac 7.540 7.490 0.050 
213 Ac 7.498 7.491 0.007 
214 Ac 7.352 7.531 -0.179 
215 Ac 7.746 7.718 0.028 
216 Ac 9.241 9.012 0.229 
217 Ac 9.832 9.931 -0.099 
218 Ac 9.384 9.437 -0.053 
219 Ac 8.826 8.818 0.008 
220 Ac 8.348 8.324 0.024 
221 Ac 7.791 7.741 0.050 
222 Ac 7.137 7.226 -0.089 
223 Ас 6.783 6.761 0.022 
224 Ac 6.327 6.318 0.009 
225 Ac 5.935 5.924 0.011 
226 Ас 5.506 5.483 0.023 
227 Ас 5.042 5.115 -0.073 
228 Ac 4.721 4.697 0.024 
229 Ac 4.444 4.382 0.062 
230 Ac 3.893 3.934 -0.041 
231 AC 3.655 3.679 -0.024 
232 Ac 3.345 3.345 0.000 
233 Ac 3.215 3.197 0.018 
234 Ac 2.930 2.942 -0.012 
235 Ac 2.852 2.886 -0.034 
236 Ac 2.723 2.668 0.055 
207Th [10] 8.328 8.277 0.051 
208TH 8.202 8.210 -0.008 
210 Th 8.069 8.065 0.004 
20TH 7.937 7.947 -0.010 
222TH 7.958 7.927 0.031 
19TH 7.837 7.817 0.020 
San 7.827 7.813 0.014 
219 hy 7.665 7.840 -0.175 
Sek? 8.072 8.056 0.016 
217Th 9.435 9.184 0.251 
218TH 9.849 9.971 -0.122 
219 Th 9.506 9.531 -0.025 
2207 8.973 8.994 -0.021 
221TH 8.625 8.595 0.030 
222 Dh 8.133 8.084 0.049 
225Тң 7.567 7.656 -0.089 
224TH 7.299 7.275 0.024 
225TH 6.921 6.884 0.037 


Table 1: (continued) 


Table 1: (continued) 


Nucl. QEPE (MeV) Theo. (MeV) AQ. (MeV) Nucl. QEPE: (MeV) Theo. (MeV) AQa (MeV) 
226TH 6.453 6.491 -0.038 2191) 9,950 9.780 0.170 
27TH 6.147 6.068 0.079 2211) 9.889 9.965 -0.076 
228 Th 5.520 5.598 -0.078 2221) 9.481 9,459 0.022 
229TH 5.168 5.124 0.044 2931) 9,158 9.113 0.045 
230TH 4.770 4.758 0.012 2247 8.628 8.580 0.048 
231TH 4.213 4.289 -0.076 22517 8.007 8.107 -0.100 
232TH 4.082 4.052 0.030 2261) 7.701 7.662 0.039 
233TH 3.745 3.757 -0.012 2271) 7.235 7.230 0.005 
234TH 3.672 3.643 0.029 2281) 6.800 6.828 -0.028 
235TH 3.376 3.406 -0.030 22907 6.476 6.413 0.063 
238TH 3.333 3.344 -0.011 2300 5.992 6.030 -0.038 
237 Th 3.196 3.146 0.050 2811) 5.576 5.608 -0.032 
211 pa 8.481 8.467 0.014 2320 5.414 5.345 0.069 
212ра 8.411 8.418 -0.007 2330 4.909 4.994 -0.085 
213pg 8.384 8.354 0.030 2310 4.858 4.860 -0.002 
214ра 8.271 8.265 0.006 29577 4.678 4.629 0.049 
215pg 8.236 8.212 0.024 2361) 4.573 4.551 0.022 
216ра 8.099 8.269 -0.170 2370 4.234 4.290 -0.056 
217ра 8.489 8.492 -0.003 23817 4.270 4.273 -0.003 
218pg 9.791 9.533 0.258 2390 4.130 4.078 0.052 
219pg 10.128 10.233 -0.105 2400] 4.035 4.067 -0.032 
220ра 9.704 9.762 -0.058 219Np 9.207 9.238 -0.031 
221 Pa 9.248 9.225 0.023 220Np 10.226 10.100 0.126 
222 Pa 8.789 8.784 0.005 222Np 10.200 10.222 -0.022 
223 Pa 8.343 8.270 0.073 223Np 9.650 9.664 -0.014 
224ра 7.694 7.788 -0.094 224Np 9.329 9.323 0.006 
225 Pa 7.401 7.379 0.022 225Np 8.818 8.765 0.053 
226 Pa 6.987 6.965 0.022 226Np 8.328 8.363 -0.035 
227ра 6.580 6.610 -0.030 227 Мр 7.816 7.847 -0.031 
228pg 6.265 6.226 0.039 229Np 7.020 7.061 -0.041 
229 pg 5.835 5.866 -0.031 230Np 6.778 6.757 0.021 
280ра 5.439 5.432 0.007 231Np 6.368 6.338 0.030 
231 pa 5.150 5.102 0.048 233Np 5.627 5.645 -0.018 
232pg 4.627 4.658 -0.031 234Np 5.356 5.376 -0.020 
233 pq 4.375 4.403 -0.028 235Npp 5.194 5.184 0.010 
234pg 4.076 4.110 -0.034 236Npp 5.007 5.021 -0.014 
235pg 4.101 4.035 0.066 237Np 4.957 4.908 0.049 
236 Pa 3.755 3.810 -0.055 238Np 4.691 4.723 -0.032 
237 Da 3.795 3.795 0.000 239Np 4.597 4.640 -0.043 
238 Pa 3.628 3.573 0.055 240Np 4.557 4.474 0.083 
2151) 8.588 8.569 0.019 24 ND 4.363 4.363 0.000 
2161) 8.531 8.570 -0.039 242Np 4.098 4.123 -0.025 
21817 8.775 8.840 -0.065 228ру 7.940 7.910 0.030 


Table 1: (continued) 


Table 1: (continued) 


Nucl. О? (MeV) Theo. (MeV) AQ. (MeV) Nucl. QEPE (MeV) Theo. (MeV) AQ. (MeV) 
229рц 7.598 7.532 0.066 "Hm 5.170 5.155 0.015 
230 py 7.178 7.207 -0.029 234 Bk 8.099 7.882 0.217 
231 py 6.839 6.890 -0.051 243 Bk 6.874 6.909 -0.035 
232 py 6.716 6.689 0.027 244 Bk 6.779 6.724 0.055 
233 py 6.416 6.426 -0.010 245 Bk 6.455 6.419 0.036 
234 Dy 6.310 6.261 0.049 246 Bk 6.074 6.149 -0.075 
235 py 5.951 6.011 -0.060 247 Bk 5.890 5.896 -0.006 
236 Pu 5.867 5.883 -0.016 218 Вк 5.827 5.765 0.062 
287рц 5.748 5.697 0.051 249Bk 5.521 5.610 -0.089 
238 Py 5.593 5.555 0.038 237 OF 8.220 8.249 -0.029 
239 Dy 5.245 5.332 -0.087 "cr 7.763 7.886 -0.123 
240рц 5.256 5.248 0.008 240CF 7711 7.745 -0.034 
241 py 5.140 5.094 0.046 кеге 7.517 7.541 -0.024 
242рц 4.984 4.982 0.002 2с 7.329 7.337 -0.008 
243рц 4.757 4.787 -0.030 245 Cf 7.258 7.169 0.089 
244 Dy 4.666 4.661 0.005 246 OF 6.862 6.862 0.000 
229 A 8.137 8.123 0.014 247 OF 6.503 6.585 -0.082 
235 Am 6.576 6.622 -0.046 248 Of 6.361 6.358 0.003 
236 Am 6.256 6.378 -0.122 2190 6.293 6.263 0.030 
238 Am 6.042 6.038 0.004 250CF 6.129 6.174 -0.045 
239 Am 5.922 5.909 0.013 251Cf 6.177 6.175 0.002 
240 Am 5.707 5.731 -0.024 292 0f 6.217 6.166 0.051 
24 Am 5.638 5.667 -0.029 253 Of 6.126 6.166 -0.040 
242 Am 5.589 5.519 0.070 254 Qf 5.927 5.915 0.012 
243 Am 5.439 5.413 0.026 241; 8.259 8.336 -0.077 
244 Am 5.138 5.207 -0.069 242 gg 8.160 8.062 0.098 
245 Am 5.160 5.152 0.008 243; 8.072 7.905 0.167 
233 Cm 7.473 7.518 -0.045 245; 7.909 7.610 0.299 
234Cm 7.365 7.382 -0.017 247; 7.464 7.378 0.086 
236 Cm 7.067 7.041 0.026 251; 6.597 6.709 -0.112 
237Cm 6.770 6.815 -0.045 252; 6.739 6.702 0.037 
238 Cm 6.670 6.676 -0.006 2536; 6.739 6.683 0.056 
239Cm 6.540 6.498 0.042 254; 6.617 6.676 -0.059 
Wée 6.398 6.396 0.002 255 Rs 6.436 6.415 0.021 
2410m 6.185 6.248 -0.063 235m 8.689 9.127 -0.438 
Mm 6.216 6.208 0.008 265m 8.379 8.391 -0.012 
M30m 6.169 6.083 0.086 M Em 8.258 8.105 0.153 
244Cm 5.902 5.910 -0.008 248 Fm 7.995 7.980 0.015 
245 Cm 5.624 5.657 -0.033 249Fm 7.709 7.713 -0.004 
246 Cm 5.475 5.489 -0.014 250 7.557 7.563 -0.006 
247 Cm 5.354 5.311 0.043 251 Bm 7.424 7.359 0.065 
248Cm 5.162 5.207 -0.045 252, 7.154 7.255 -0.101 
249Cm 5.148 5.154 -0.006 253 ту 7.198 7.192 0.006 
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Table 1: (continued) 


Nucl. О? (MeV) Theo. (MeV) AQ. (MeV) 
254 7.307 7.256 0.051 
255р 7.241 7.259 -0.018 
256 7.025 7.032. -0.007 
257 Ет 6.864 6.882 -0.018 
%6М4 8.889 9,193 -0.304 
MTM 8.764 8.983 -0.219 
зама 8.497 8.647 -0.150 
250Md 8.155 8.135 0.020 
251Md 7.963 7.982 -0.019 
253Md 7.573 7.814 -0.241 
255Md 7.906 7.834 0.072 
257Md 7.557 7.505 0.052 
258Md 7.271 7.263 0.008 
251 No 8.752 8.833 -0.081 
25230 8.549 8.555 -0.006 
25330 8.415 8.406 0.009 
25430 8.226 8.327 -0.101 
25530 8.428 8.413 0.015 
25630 8.582 8.480 0.102 
257No 8.477 8.496 -0.019 
259No 7.854 7.859 -0.005 
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Due to the successful calculations of the binding energies, 
it is expected that the a-decay energies, which are the dif- 
ferences among the binding energies of the parent nuclei, 
the daughter nuclei, and the a-particles, can be reproduced 
with good accuracy. Thus, we calculate the a-decay ener- 
gies for 1169 nuclei with 50 < Z < 110 and compare 
the calculated results with the experimental data taken from 
AME2020 [27]. The deviations between the experimental 
a-decay energies and the theoretical results for these nu- 
clei are depicted in Fig. 3. In Fig. 3, the blue circles de- 
note the deviations and the red shadow shows the deviations 
WEE - Сы < 0.3 MeV. The dashed line represents 
GEN - Ge = 0 MeV and the two dash dotted lines 
present igor E gue — 0.5 MeV, respectively. The de- 
viations for the o-decay energies of the 1169 nuclei are all 
clearly below 0.5 MeV and the deviations for most of these 
nuclei are less than 0.3 MeV. These results show good agree- 
ment between the theoretical a-decay energies derived from 
the binding energies which are calculated using the Gaussian 
process and the experimental data. Furthermore, it has been 
found in previous studies that a-decay energies are strongly 
affected by the shell effect, which leads to larger deviations 
for nuclei near the closed shell [56]. In Fig. 3, the deviations 
for nuclei near the shell closure are also less than 0.3 MeV. 
It can reflect that 7 and v features can successfully model 
the shell effect with the Gaussian process. We also calculate 
the average deviation and standard deviation for these nuclei, 
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Fig. 3. The deviations between the experimental a-decay energies 
and the theoretical results for 1169 nuclei with 50 < Z < 110. 
The blue circles depict the deviations for these nuclei. The dashed 
line denotes on — poe = 0 MeV. The red shadow and 


the dash dotted lines present [quoe - oe | < 0.3 MeV and 
los" = Qa] = 0.5 MeV, respectively. 


given by 


Па 


= ES Expt.p _ rjTheo.p| _ 
(Ta) = fig 2-10 mp — QTheo.p| = 0.047 MeV (7) 


and 
i 2 
02 = |У (qr z SEH = 0.070 MeV. 
Па 
р=1 
(8) 


Owing to the complexity of ће quantum many-body theory, 
itis difficult to calculate the a-decay energies with deviations 
less than 0.1 MeV. These small deviations show that the a- 
decay energies agree well with the experimental results. 
Recently, some actinide nuclei, including 204Ас [57] and 
207Th [10], were synthesized experimentally. Theoretical a- 
decay properties provide useful references for these experi- 
ments. Here, we present the theoretical a-decay energies cal- 
culated using the Gaussian process for the actinide nuclei in 
Table 1. In Table 1, the first column lists the actinide nu- 
clei. The second column denotes the experimental data and 
the third column presents the theoretical results. The fourth 
column gives the deviations AQ, = QE*Pt. — (Theo. pe. 
tween the experimental results and the theoretical ones. The 
experimental a-decay energies for two new nuclides 204Ас 
and 2°’ Th are taken from Ref. [57] and Ref. [10] separately. 
It can be clearly seen that the theoretical results obtained us- 
ing the Gaussian process are in good agreement with the ex- 
perimental data for the actinide nuclei. For the new nuclide 
204 Ac, the theoretical a-decay energy calculated using the 
Gaussian process is nearly equivalent to the experimental re- 
sult, with a small deviation of AQ, = (ЕР — (9Theo. = 
0.0004 MeV. For another new nuclide, 2°’Th, the deviation is 


ж Аба = QExPt. — QTheo. — 0.051 MeV, indicating that the 
calculated result is in good agreement with the experimental 
one. These results demonstrate that the a-decay energies for 
the actinide nuclei can be well reproduced by deriving from 
аз the theoretical binding energies calculated using the Gaussian 
414 process. Overall, the above results show the reliability of the 
Gaussian process in the calculations of nuclear binding ener- 
gies and a-decay properties. 
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Table 2: The predicted a-decay energies and half-lives for some un- 
known actinide nuclei. The first column denotes the a-decay emitters. 


Table 2: (continued) 


Nucl. 98 (MeV) ФАРМ (MeV) logio(T55) guef" 


The second and third columns are the predicted a-decay energies cal- 


culated using the Gaussian process and the FRDM separately. The 
fourth and fifth columns represent the predicted a-decay half-lives cal- 
culated using the NGNL with the predicted a-decay energies given by 
the Gaussian process and the FRDM, respectively. The units of the 


a-decay half-lives are seconds. 


Nucl. ОС” (MeV) ФХР“ (MeV) logio(Tp3) logio(T13 


200 Ac 9.260 8.905 -4.982 -4.089 
201 Ac 9.016 8.895 -4.373 -4.062. 
202 Ac 8.639 8.685 -3.383 -3.507 
203 Ac 8.432 8.575 -2.811 -3.208 
203TH 8.948 8.825 -3.865 -3.542 
204TH 8.827 8.765 -3.547 -3.381 
205TH 8.595 8.575 -2.917 -2.862 
206Th 8.469 8.515 -2.565 -2.694 
207ра 8.495 8.765 -2.289 -3.040 
208ра 8.478 8.565 -2.240 -2.487 
209ра 8.461 8.305 -2.191 -1.738 
210ра 8.456 8.265 -2.176 -1.619 
21917 8.456 8.605 -1.825 -2,252 
210] 8.512 8.485 -1.986 -1.909 
arei 8.490 8.365 -1.922 -1.558 
21911 8.542 8.385 -2.071 -1.617 
215Np 8.444 8.815 -1.435 -2.490 
216Npp 8.544 8.625 -1.726 -1.958 
217Np 8.684 8.725 -2.124 -2.239 
218Np 8.956 8.945 -2.872 -2.842 
224Pu 9.914 9.565 -5.944 -5.107 
225 pg 9.306 9.285 -4.455 -4.401 
226 Pu 8.774 9.035 -3.028 -3.744 
227рц 8.246 8.695 -1.476 -2.804 
225 Am 9.977 9.895 -5.779 -5.587 
226 Am 9.348 9.605 -4.235 -4.884 
227 Am 8.936 9.345 -3.136 -4.227 
228 Am 8.423 9.075 -1.657 -3.515 
290m 8.727 9.395 -2.202 -4.027 
230Cm 8.419 8.725 -1.288 -2.196 
231 Cm 8.068 8.385 -0.182 -1.183 


232 Cm 7.788 7.885 0.753 0.424 
229Bk 9.233 9.555 -3.270 -4.112 
230 Bik 8.864 9.165 -2.249 -3.086 
231Bk 8.549 8.805 -1.326 -2.080 
232Bk 8.270 8.465 -0.464 -1.070 
2380 9.392 8.585 -3.360 -1.080 
2M E 9.092 8.665 -2.548 -1.319 
295 OF 8.726 8.535 -1.500 -0.927 
2360 8.509 8.335 -0.847 -0.305 
237; 9.679 8.645 -3,778 -0.906 
238; 9.215 8.485 -2.548 -0.415 
) 239Rs 8.935 8.155 -1.760 0.643 
240 gs 8.558 7.975 -0.639 1.248 
239 Eg 10.540 8.845 -5.536 -1.152 
0 10.260 8.605 -4.887 -0.428 
24 АЕ 9.791 8.405 -3.739 0.199 
242Fm 9.548 8.285 -3.110 0.586 
M2Mg 10.392 9.045 -4.889 -1.389 
M3Mg 10.139 9.005 -4.284 -1.273 
“ама 9.767 8.935 -3.854 -1.068 
%5ма 9.541 8.925 -2.762 -1.038 
245No 10.326 9.505 -4.423 -2.336 
246No 10.104 9.465 -3.883 -2.227 
247No 9.878 9.335 -3.316 -1.869 
248No 9.629 9.205 -2.667 -1.504 


Finally, we predict the a-decay energies for some unknown 
ав actinide nuclei using the Gaussian process. With the pre- 
4:9 dicted a-decay energies, we also calculate the a-decay half- 
42 lives using the new Geiger-Nuttall law (NGNL) [58]. The 
41 Corresponding results are given in Table 2. In Table 2, the first 
42 column lists the a-emitters. The second and third columns 
45 present the a-decay energies calculated using the Gaussian 
44 process and the FRDM, respectively. The fourth and fifth 
45 columns give the predictive o-decay half-lives calculated us- 
а ing the NGNL with the o-decay energies predicted by the 
47 Gaussian process and the FRDM, respectively. It can be 
ав found that most predicted a-decay energies agree well with 
429 those calculated using the FRDM. Nevertheless, the predicted 
азо 0-ӣесау energies for Einsteinium, Fermium, Mendelevium, 
431 and Nobelium are relatively larger than those given by the 
42 FRDM, which results in different a-decay half-lives. We 
433 hope that future experimental a-decay properties for Ein- 
аза Steinium, Fermium, Mendelevium, and Nobelium can pro- 
аз Vide useful information for improving the Gaussian process. 
аз The a-decay properties predicted by the Gaussian process can 
437 complement existing theoretical models and provide valuable 
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tinide isotopes are being synthesized at the Heavy Ion Re- 
search Facility in Lanzhou (HIRFL), China. Therefore, it is 
expected that the predicted a-decay properties can be used 
as theoretical references for identifying new nuclides in the 
future. 


IV. SUMMARY 


In this work, the Gaussian process with a composite kernel 
function is applied to study the binding energies. First, we 
calculate the binding energies for 2238 nuclei with Z > 20 
and N > 20 within the framework of the Gaussian process 
using a physically motivated feature space. The calculated 
average deviation and standard deviation are 0.046 MeV and 
0.066 MeV, respectively. The results demonstrate that the 
binding energies are successfully modeled by the Gaussian 
process, reflecting the good learning ability of the Gaussian 
process in the calculations of binding energies. Then, we cal- 
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